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ABSTRACT 

Secondary contributions to the anisotropy of the Cosmic Microwave Background (CMB), such as the 
integrated Sachs-Wolfe (ISW) effect, the thermal Sunyaev-Zel'dovich effect (tSZ), and the effect of 
gravitational lensing, have distinctive non-Gaussian signatures, and full descriptions therefore require 
information beyond that contained in their power spectra. The Minkowski Functionals (MF) are well- 
known as tools for quantifying any departure from Gaussianity and are affected by noise and other 
sources of confusion in a different way from the usual methods based on higher-order moments or 
polyspectra, thus providing complementary tools for CMB analysis and cross-validation of results. In 
this paper we use the recently introduced skew-spectra associated with the MFs to probe the topology 
of CMB maps to probe the secondary non-Gaussianity as a function of beam-smoothing in order to 
separate various contributions. We devise estimators for these spectra in the presence of a realistic 
observational masks and present expressions for their covariance as a function of instrumental noise. 
Specific results are derived for the mixed ISW-lensing and tSZ-lensing bispectra as well as contam- 
ination due to point sources for noise levels that correspond to the Planck (143 GHz channel) and 
EPIC (150 GHz channel) experiments. The cumulative signal to noise ration S/N for one-point gen- 
eralized skewness-parameters can reach an order of O(10) for Planck and two orders of magnitude 
higher for EPIC, i.e. 0(1O 3 ). We also find that these three spectra skew-spectra are correlated, having 
correlation coefficients r ~ 0.5 — 1.0; higher I modes are more strongly correlated. Though the values 
of S/N increase with decreasing noise, the triplets of skew-spectra that determine the MFs become 
more correlated; the S/N ratios of lensing-induced skew-spectra are smaller compared to that of a 
frequency-cleaned tSZ map. 



1 INTRODUCTION 



All-sky multi-frequency Cosmic Microwave Background (C MB) missions, such as the completed WMAlffl ongoing Plancljf] and future (proposed) 



Experimental Probe of Inflation ary Cosmology (EPIC) survey teock et alj2008 . 2009; B aumann et al.l2009l) or ESAs Cosmic Origin Explorer (COrE, 
iThe COrE Collaboration! J2012h ) are major sources of information about the properties of the primordial density fluct uations that seeded the process 
of ga l axy formation in the Un iverse as well as other key asp ects of cosmological theory, including the g lobal isotropy dCopi et al]|2007l ; lHoftuft et al.l 
l2009l ; lHanson & Lewisll2009h and topology of the Universe l lLuminet et al]2003l ; lRoukemia et alj2004l) . 

The study of non-Gaussianity in the CMB fluctuations can provide valuable and detailed information regarding the physics of the early Universe 
of the inflationary epoch In the standard slow-roll paradigm, the scalar field responsible for inflation fluctuates with a minim al amount of self 



interaction which ensures tha t any non-Gaussianity generated during the i nflation through s elf-in teraction is expected to be sm all I Salopek & Bond 
Il990l. Il99ll ; iFalk et al.l [l993l : iGangui et alj|l994l ; kcquaviva et~al]|2003l ; lMaldacenj|2003l) : see iBartolo. Matarrese & Riottcl J2006l) for a review. 
Variants of the simple inflationary model such as multiple scalar fields, features in the inflationary potential, non-adiabatic fluctuations, non-sta ndard 
kineti c terms, warm inflation, or deviations from Bunch-Davies vacuum can however all lead to higher level of primordial non-Gaussianity dCherJ 
l2010h . 

However, the detection of departure from Gaussianity in the CMB can be due to either primary or secondary effects (or both), as well as the mode- 
coupling effects of secondaries and gravitational lensing along the observer's light cone. Secondary anisotropics resulting from the formation of struc- 

1 http://map.gsfc.nasa.gov/ 

2 http://www.rssd. esa.int/index.php?project=Planck 
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ture are known to dominate at smaller angular scales, are highly non-Gaussian in nature dCooravll20oibl ; ICoorav & Hull200d : lve"rde & Spergelll2002l) 
and are arguably as interesting as their primary counterpart. One of the prominent contributions to the secondary non-Gaussianit y is due to the mode- 
coupli ng of weak gravitation al lensing and sources of secondary contributions such as the thermal Sunyaev-Zel'dovich effect dGoldberg & SpergeJ 
ll999d ltt lCoorav & Hvj|200ch . Although weak lensing of the CMB produces its own characteristic signature in the angular power spectrum, its de- 
tection has proved to be difficult using the CMB power spectrum alone. Non-Gaussianity imprinted by lensing into the primordial CMB remains 
below the detection level of current experiments, although with Planck the situation is likely to improve. Nevertheless, cross-correlating CMB data 
with external tracers means lensing signa ls can be probed at th e level of the mixed bispectrum. After the first unsuccessful attempt to cross-correlate 
WMAP against SDSS, recent efforts by ISmith. Zahn & Dord d2007t) have found a clear signal of weak lensing of the CMB, by cross-correlating 
WMAP against NVSS. Their work also underlines the link between three-point statistical estimators and the estimators for weak lensing effects on 
CMB. The understanding of secondar ies are not only important in their own right, but also from the perspective of their impact on estimation of 

cosmological parameters dSmidfeOld) . 

The study of non-Gaussianity is usually primarily focused on the bispectrum, as this saturates the Cramer-Rao bound dBabichl 120051 : 
iKamionkowski. Smith & Heavens . I [20 llh and is therefore in a sense optimal, ho wever in practice it is diffi cult to probe the entire configuration 
dependence in harmonic space contained within the bispectrum using noisy data dMunshi & HeavensllioTol) . The cumulant correlators are multi- 
point correlators collapsed to encode two-point statistics. These were introduced in the context of analyzing g alaxy clustering bvlSzapudi & Szalavl 
and were later found to be useful for analyzing projected surveys such as the APM galaxy survev dMunshi. Melott & Coles|l200fj) . Be- 
ing two-point statistics they ca n be analyzed in multipole space by defining an associated power-spectrum. Recent studies bv lCoorav and 
ICoorav. Li & Melchiorril d2008h have demonstrated their wider applicability including, e.g ., in 21cm studies. In more re cent studies the skew - and 
kurt-spectra were found to be useful for analys ing temperatu re I Munshi & Heavensll2010l) as well as polarization maps dMunshi etai]|201 id) from 
CMB experiments and in weak lensing studies dMunshi et alfeOllbl ld). 

In addition to studies involving lower order multi-spectra, MFs have been extensively developed as a statistical tool for non-Gaussianity in a cos- 
mological setting for both 2-dimensional (projected) and 3-dimensional (redshift) surveys. Analytic results are known certain properties of the MFs 
of a Gaussian random field making them suitable for identifying non-Gaussianity. Exam ples of such stu dies include CMB data |Schmalzing & Gorskil 



1998l:lNovikoy. Schrnalzing and Mukhanoyll2o6ol ; lHikage et alj2008tlNatoli et alj201Ch. weak lens i ng dMatsubara and Jainld200ll);ISato et al. 
Taruva et al.l d2002h: IMunshi et al.l d2012ah). large-scale structure dGott et alj|l986l: IColesI 1 19881: IGott et alj|l989l: lMelotdll99cl: iGott et al 



20011) ; 



19901: 



Moore et aUll992t Gott et alll 19921: JCanavezes et alj|l998l; ISahni, Sathvaprakash & Shandarin|ll998l ; ISchmalzing & Diaferioll2000l : iKerscher et all 



Gott et al.L 

20011 ; iHikage et alJl2002L IPark et alj|2005l: IHikage. Komatsu & Mastubarall2006l : lHikage et alj|2008h. 21cm dGleser et alJl200d)Tlre quency cleaned 
Sunyaev-Zel'dovich (SZ) maps dMunshi et all 1201 2bl) and N-body simulations dSchmalzing & Diaferiol 120001 : IKerscher et alj|200lh . The MFs are 
spatially-defined topological statistics and, by definition, contain statistical information of all orders in the moments. This makes them complemen- 
tary to the poly-spectra methods that are defined in Fourier space. It is also possible that the two approaches will be sensitive to different aspects of 
non-Gaussi anity and systematic effects although in the weakly non-Gaussian limit it has been shown that the MFs reduce to a weighted probe of the 
bispectrum dHikage. Komatsu & Mastubar j f2006l) ). 

The skew-spectrum is a weighted statistic that can be tuned to a particular form of non-Gaussianity, such as that which may arise either during 
inflation at an early stage or from structure formation at a later time. The skew-spectrum retains more information about the specific form of non- 
Gaussianity than the (one-point) skewness parameter alone. This allows not only the exploration of primary and secondary non-Gaussianity but also 
the residuals from galactic foreground and unresolved point sources. The s kew-spectrum is directly related to the lowest order cumulant correlator 
and is also known as the two-to-one spectra in the literature dCooray||200T3) . In a series of recent pu blications the concept of skew-sp e ctra was gen- 
eralized to analyse the morph o logica l properties of cosmological data sets or in particular the MFs bv lMunshi. Smidt & Cooravl d20ld) ; lMunshi et al] 
d2012dlbT); IPratten & MunshU d2012h . The first of these three spectra, in the context of secondary-lensing cor relation studies, was introduced by 
IMunshi et alj 1 20 Hal) and was subsequently used to analyse data release from WMAP bv lCalabrese et all d20ld) . 

The primary aim of this paper is to consider the entire set of generalised skew-spectra resulting from the mode-coupling of secondary anisotropies 
and lensing of the CMB and the contribution thereof to non-Gaussian morphology of the CMB maps. We will be considering three different secondary- 
lensing correl ation bispectra. The secondaries that we consider are the Integrated Sachs-Wolfe effect ( ISW) that dominat es at large angular scales 
(Cooray 2002) and the thermal Sunyaev-Zel'dovich (tSZ) effect that dominates at smaller angular scales iBirkinshawl d 19991) . In addition we consider 
a foreground, namely the contribution from unresolved point sources. We will consider two experimental setups, the ongoing Planck satellite and the 
the proposed EPIC satellite mission discussed above. 

The layout of the paper is as follows. In Sj2]we briefly outline the bispectrum corresponding to lensing-secondary mode-coupling. Next, in ij3] 
we review the formalism underlying the Minkowski Functionals and in Sj4]we introduce the generalised skew-spectra associated with the MFs. In Sj5] 
we present the estimators for these spectra and their covariances. Fin ally, in £]7]we discus s our results and comment on future implementation. 

Throughout we will use the parameters of a WMAP cosmology dLarson et alj|201 ll) . 



2 MODE COUPLING INDUCED BY LENSING - SECONDARY CROSS-CORRELATION AND THE RESULTING BISPECTRUM 

The bispectrum of primary anistropies encodes information that can be used to constrain the inflationary dynamics but, as discussed in the previous 
section, the primary contribution to non-Gaussianity is expected to be negligible in the simplest realisations of the generic inflationary scenario. 

The secondary bispectrum provides valuable information regarding the low-redshift Universe and constrains structure formation scenarios. The 
secondaries can be broadly divided into three different types: 
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Figure 1. The cross-spectra ft for various secondaries and lensing are plotted as a function of the harmonics I. From left to right different panels correspond to cross- 
correlation of lensing potential and SZ, ISW and PS contribution. The cross-spectra /3i is being used in Eq.(5J for the construction of mixed bispectrum -E^^ 2 S ; 3 ■ It is 
defined in Eq.jTJ. Various estimators for the skew-spectra that we will use, associated with the three MFs, will be define d using the mix ed bispectra. A background 
ACDM cosmology is assumed. The details of these calculations, which relies on halo model prescriptions, can be found in iCooravlboOlal) . 



(i) Gravitational secondaries, ca used by evolution in the gravitational potential along the observer's past light cone including the well-known 
integrated Sachs-Wolfe (ISW) effect dKofman & Starobinskvll985l:|Martinez-Gonzalez. Sanz. Silkll99d:lMukhanov. Feldman & BandenbergeJ 19921 : 
iKaminkowski & Spergellll994 iMunshi. Souradeep & Starobinskvlll995l : iMollerach et alj|l995t Boughn & Crittendenll2004l) as well as the Rees- 
Sciama (RS) effect. 

(ii) Scattering secondaries, such as the thermal Sunyae v-Zerdovich (tSZ) effect teirkinshawll 19991) . kinetic Sunyaev Zel'dovich (kSZ) effect and 
the Ostriker-Vishniac effect (see e.g. ICastro et al.l d2004l )). These effects are caused by the interaction of the CMB photons with the intervening 
free-electron population. 

(iii) Lensing secondaries caused by the propagation of photons through large scale structures. 

Contributions to secondary bispectra can also arise from terms involving the cross-correlation of gravitational lensing and the effects of inter- 
vening material, such as the tSZ effect due to inverse Compton scattering of CMB photons from hot gas in the intervening clusters. The decay of 
the peculiar gravitational potential along the line of sight in ACDM cosmology, introduced above as the Integrated Sachs- Wolfe or ISW effect, is 
correlated to the lensi ng due to the potential, can also generate an additional contribution to the secondary bispectrum in a similar fashion; see e.g. 
ICoorav & Sett] d2002h for a detailed discussion of various secondaries in the context of halo model. The contribution to secondaries due to reioniza- 
tion of the Universe are detailed in lHu. Scott & Silkl d 19941) . Foregrounds, such as unresolved point sources (PS), can also contribute to the secondary 
bispectrum through their cross-correlation with the lensing of CMB. 

On a different note, we comment that while the study of secondary anisotropies i s important in their own right, they are also important in 
their effect on the calculation of error covariances in cosmological parameter estimation JjoudakietalfeoiOh . Understanding the detailed statistical 
properties of secondary anisotropies like those we discuss here is therefore is of the utmost importance in the era of precision cosmology. 

We will be dealing w i th the secondary bi s pectra invol ving the lensing of both primary anistropies and other secondaries. Following 
iGoldberg & Spergell Jl999ah . iGoldberg & Spergell Il999bl) and ICoorav &~Hv] bOOfj) we start by expanding the observed temperature anisotropy 
in terms of the primary contribution Op(Cl), the secondary contribution Bs(O) and lensing of the primary Ol(Q,): 

e(fi) = e P (n) + e L (A) + e s («) + ••• . (i) 

Here Cl — (8, <f>) is the angular position on the surface of the sky. Expanding the respective contribution in terms of spherical harmonics Yi m (Cl) we 
can write: 

G P (n) = ^(e p ) im y !ro (f2); G L (n)s^[w(fi)-ve P (n)] !m Fi m (n); e s (n) = ^2(@ s )imYi m (Ci). (2) 

lm Im Ira 

Here ip(Cl) is the projected lensing potential ( IGoldberg & Spergell 1999al lbl). The second ary bispectrum for th e CMB takes contributions from products 
of P, L and S terms with varying order. The bispectrum Bf^, are defined as follows (see lBartolo et al.l d2004l> for generic discussion on the bispectrum 
and its symmetry properties): 

B lhh = E ( I 1 I 2 t ) /(0p(ni)0L(«2)e s (n3))y ! : mi (ni)yi; ro2 (fi2)y ! ; m3 (n3)dnidn 2 dn 3 ; 

V TTb \ TTL'2, Tib ,3 J J ' 

((Sp)i 1 m 1 (SL)( 2 m 2 (Ss)l 3 m 3 )- (3) 



m\ m<2 m3 



mi ni2 m3 
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The angular brackets represent ensemble averages. The matrices denote 3j symbols ( lEdmonds|[l963) and the asterisks denote complex conjugation. 
It is possible to invert the relation assuming isotropy of the background Universe: 

<(ep) iimi (e L ) i2m2 (e s ) i3m ,) = ( ^ ^ l ^) B ™ 3 - W 

Finally the bispectrum Bf^i is expressed in terms of the un-lensed primary power spectrum Ci = ((0p)i m (6p)!m) and the cross-spectra /3i (to be 
defined below) as follows: 



B hl2h = ~ { Ph C h ' ^' " 1( " ± ri ' " JV " J ' ^' + eye. perm. J- h t i 2 i 3 = Bi^iji^i^ 



Uh + 1) - h(h + 1) - h(h + 1) 



(2h + 1)(2Z 2 + l)(2l 3 + 1) / h h h 



(6) 



4tt V 

(see lGoldberg & Spergell Jl9993) . I Goldberg & Spergell ll999bh for a derivation). The reduced bispectrum above is denoted Bi x i 2 i 3 . To simplify the 
notation for the rest of this paper, we henceforth drop the superscript PLS form the bispectrum Bi 1 i 2 i 3 . The cross-spectrum /3; introduced above 
represents the cross-correlation between the projected lensing potential tp(Cl) and the secondary contribution 0s(O): 

Mft)e s (n')> = — £)(2j + i)#Pt(fi.n'). (7) 

1=2 

The cross-spectra /?; take different forms for ISW-lensing, RS-lensing or SZ-lensing correlation and we assume a zero primordial non-Gaussianity. 
The reduced bispectrum Bi t i 2 i 3 defined above using the notation /j^i, is useful in separating the angular dependence from the dependence on power 
spectra C; and /?;. We will use this to express the topological properties of the CMB maps. The parameters for lensing secondary correlations are 
displayed in Figure-Q] The left, middle and right panels in Figure-CD displays SZ-lensing, IS W-lensing and point source lensing correlations. These 
results are based on halo model calculations performed using the halo model l lCooravll2001al) . 

The beam bi (9b) and the noise of a specific experiment are characterised by the parameters erbeam and <r r 



' rms ■ 



bl(9 b ) = exp[-l(l + l)Obeam]; GTbeam = b ; 71; = ^ mB Qpi x \ fipix = r^- . (8) 

i/81n(2j i^pix 

where a Tms is the rms noise per pixel that depends on the full width at half maxima or FWHM of the beam 9b- The number of pixel N p i x required 
to cover the sky determines the size of the pixels f2 p i x . To incorporate the effect of experimental noise and the beam we have to replace Ci — > 
Cibf(6b) + ni, and the normalization of the skew-spectra that we will introduce later will be affected by the experimental beam and noise. The 
computation of scatter will also depend on these parameters. We will consider two different experimental setups: Planck and EPIC. The parameters 
of these experiments are tabulated in Table[Q 

The optimal estimators for lensing-secondary mode-coupling bispectrum have been recently discussed bv lMunshi et alj l feoilal) . The estimators 
that we propose here are relevant in the context of constructing the MFs. 



3 MINKOWSKI FUNCTIONALS 

Integral geometry provides a natural framework within which to define the set of morphological descriptors for a random field. These descriptors 
are intrinsically defined in the spatial domain where they take into account all n-point correlators up to arbitrary order. Hadwiger's characterization 
theorem shows that a li near comb i nation of these d + 1 functionals will provide a complete morphological description of the morphology of d- 
dimensional objects; see lHadwigej ( fl95^) for a formal treatment. These functionals are more commonly referred to as the Minkowski functionals. The 
Minkowski Functionals are usually calculated using vo lume-weighted curvature integrals for which the analytical results for a Gaussian random field 
are known (Adler 1981; Tomita 1986; Gott e t aljl990l) . More recently the analytical values for wea kly non-Gaussian fields have been calculated as a 
funct i on of skewness parameters by using a perturbative appro ach based on the Edgeworth expansion l lMatsubaral 1 99 4 , 1995; Matsubara & Yokohamal 
Il996l ; lMatsubarall2002l : iHikage. Komatsu & Mastubarall2006t) ). This approach allows us to use the MFs as a test of non-Gaussianity in the weakly 
perturbed regime as constrained by observation and predicted by models for inflation. 

In 2 dimensions the MFs Vo(u), Vi(y) and V2{y) correspond respectively to the area of a set E, length of the perimeter of the set and the 
integrated curvature along its boundary. The MF V^iy) can be related to the well-known genus g and the Euler characteristic \: 

Vo(u) = f da; Vi(u) = -f dl; V 2 {v) = ^- [ ndl. (9) 

Here dl and da represent the length and surface element respectively. In our analysis we consider a smoothed random field Q(Q) with mean (G(O)) = 
and variance cro(#&) — (@ 2 (i^)). For a generic 2-dimensional weakly non-Gaussian random field on the surface of the sky, the spherical 
harmonic decomposition using YJ m (f2) as basis functions 0(A) = QimYi m (Cl) can be used to define the power spectrum Ci which is sufficient 
to characterize an isotropic Gaussian field (0; m 6*, m ,) = CiS U iS mm i . 
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Figure 2. The variances cr (9 b ) ar| d 01 ($()), defined in Eq. )12t . are plotted as a function of the FWHM 9 b . The left panel corresponds to an experimental setup such as 
EPIC and the right panel corresponds to Planck-type experiment. See Table{T]for detail specifications regarding the noise level and beam. 



We will b e studying the MFs defi ned over the surface of the celestial sphere but equivalent results can be obtained in 3D u sing a Fourie r 
decomposition jPratten & Munshill2012l) . The MFs for a 3D random Gaussian field are well known and are given by Tomita's formula dTomital 19861) . 

For a non-Gaussian field the higher order statistics such as bi- or tri-spectrum can describe the resulting mode-mode coupling. Alternatively 
topological measures such as the MFs (including the Euler characteristic or genus) can be employed to quantify deviations from Gaussianity. Indeed 
it can be shown that the information content in both descriptions is equivalent in that, at leading order, the MFs can be constructed completely from 
the knowledge of the bispectrum alone. 

The notations and analytical results in this section are being kept generic however they will be specialized to the case of CMB sky in subsequent 
discussions. 



The MFs denoted as Vk(y) for a threshold v = O/cto; where u (8 b ) = (0 ) are perturbatively expressed 



as: 



V k (y) = (27r)( l +1)/a ^ exp (-£) (-^-) ' [vW(»)a oi e b) + Vf\vy o (0 b ) + V<»*) + •• •] (10) 

vfV) = 



\s i0 \e b )H k+2 {u) + ^s {1 \e b )H k {v) + fc(fc 6 l) s {2 \e b )H k - 2 {v) 



(ii) 



= 7- Y,( 21 + !)[*(' + l)] J Cibf(9 b ) (12) 

( 

The constant cu k introduced above is the volume of the unit sphere in k dimensions. w k = ir k ^ 2 /T{k/2 + 1) in 2-dimension we will only need 
u>o = 1, U)i = 2 and u 2 = it- Here T is the the gamma function. The lowest-order Hermite polynomials Hk(y) are listed below. 

H-i{v) = y|exp erfc (^=) ; H (v) = l, H^u) = u, 

H 2 {u) = v 2 -l, H 3 {v) = u 3 - 3u, H A {v) = v 4 - 6j/ 2 +3, 

^M = (-irexp(^)-Ae X p(4). (13) 

The expression consists of two distinct contributions. The part that does not depend on the three different skewness parameters 
S^ t) \8 b ),S ( ' 1 \9 b ),S i - 2 \8 b ) and signifies the MFs for a Gaussian random field. The other contribution 8Vk(y) represents the departure from the 
Gaussian statistics and depends on the generalized skewness parameters defined in Eq.dl5t - Eq,dl7t. Various second-order moments Uj (9b) defined 
in Eq,dl2t appear in Eq.dlOb and Eq.dl It can be expressed in terms of the power spectra Ci and the observational beam bi(8 b ), assumed Gaussian 
with a full width at half maximum 8 b \ see Eq.([8} for definition. The moment ao(8 b ) is a special case which relates to the variance of the field. The 
quantities ai(8 b ), a 2 (8 b ) are natural generalization of this variance, putting greater weight on higher-order harmonics; the variances that appear most 
frequently henceforth are al{8 b ) = (8 2 ) and cr 2 ^,,) = ((V9) 2 ). 

Real space expressions for the triplets of skewness (8 b ) are given below. These are natural generalizations of the ordinary skewness (8 b ) 
that is used in many cosmological studies. They are all cubic statistics but are constructed from different cubic combinations. 

_ <f 3 > ■ 5 U^= 3 g (92v2e) (fl t ) 3 (e 2 V 2 9) , 



%) a 4 (9 b )' 4 al{8 b )al{8 b ) 4a 2 (fl 4 KW 

aw,^- o g (ve ' v9v2e) (^) o ((ve).(ve)(v 2 0)) 

s {eb) = - 3 — — =^ 3 — — • (14) 
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Figure 3. Same as previous plot but for Planck type noise. The normalisation of the skew spectra is depends on the noise level. The shape of the bispectrum is independent 
of the noise but depends on the experimental beam. Notice that S"(9j) changes sign near 9 b ~ 10' as before. is negative for all three effects. The other 

parameters remain positive for the entire range of FWHM probed. 
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Figure 4. The absolute magnitudes of the one-point skewness parameters S^°'(9 b ), S^(9 b ) and S^(9 b ) (solid, short-dashed and long dashed lines respectively), 
are plotted as a function of smoothing scale 9 b . These parameters are defined in Eq. )14t . Various panels correspond to contributions from different types of secondary 
non-Gaussianity that represent the cross-correlation of lensing and Sunyaev-Zel'dovich effect (left-panel), cross-correlation of ISW and lensing (middle-panel) and 
cross-correlation of point sources and lensing (right-panel). The parameter S^- ' (6),) is negative for the entire range of FWHM or 9 b considered, for all three bispectra 
we have studied. The parameter S'^' (9 b ) shows a transition from being negative to positive for the case of ISW bispetcrum. All other skewness parameters stays positive 
as indicated. An experimental set up which is same as EPIC was considered. See Tablef_T]fo r detailed specifications regarding level of noise and beam. The bispectrum 
used in our calculation is given in Eq.^3) and the cross-spectra are plotted in FigureJJ] 



Notice that a knowledge of the S^(6 b ) parameters completely specifies the V fe ' 2 ' ' (v) parameters which characterize the lowest order departure 
from Gaussianity. The expressions in the harmonic domain are more useful in the context of CMB studies where these skewness parameters can 
be recovered from a masked sky using analytical tools that are commonly used for power spectrum analysis. The skewness parameter 5 , ' 1 '(0j) is 
constructed from the product field [0 2 ] and [V 2 0], whereas skewness parameter S^ 2 \6 b ) relies on the combination of [VQ ■ V0] and [V 2 0]. By 
construction, the skewness parameter S 1 ' 2 ' (6b) has the highest weight for high I modes and s£ ' has the lowest weight from high / modes. 
The expressions in terms of the bispectrum Bi 1 i 2 i 3 (see Eq.l[3} for definition) take the following form: 

s l@3) (e b ) = -]T B lll2l3 i llhl3 b h (e b )b l2 (e b )h 3 m, (i5) 

(«=2 

s (e2v2e) (<9 6 ) = __ 2 [zi (Zi + 1) + h(h + 1) + h(h + i)] B llhh i lll2l3 b h (e b )b l2 (e b )b h (e b ), (16) 

1 Imam 

^e-vev'e)^ = g + ^ + ^ + ^ + + ^ + C y C p erm J B^i^I^i^b^ (6 b )h 2 (Ob)h 3 (9b)- (17) 

i i= 2 
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Figure 5. The three skew-spectra corresponding to the three MFs (Of,) (left-panel) , Sj L> (9b) (middel-panel) and S^ 1 (9b) (right-panel) are plotted for the mixed 
secondary bispectrum of SZxLensing as a function of the harmonics I. A background ACDM cosmology is assumed. The skew-spectra are defined in Eq(BJ, Eqfi6t 
and Eq{T7) respectively. All other sources of non-Gaussianity are ignored. Four different Gaussian-beams are considered. The curves from left to right correspond to 
Qb = 10', 25', 50', 100' in each panel. The normalization of the skew spectra is somewhat arbitrary and do not affect the signal to noise ratios. We have ignored the 
presence of noise, as defined in Eq. U2t . in our calculation of co (8b) and a\ (9b) respectively. Notice that for a given 9b higher order skew-spectra peak at a higher I. 
The signal to noise of the skew-spectra associated with SZ effect are plotted in Figure|9]for EPIC and Figure[8]for Planck noise level respectively. The skew-spectra are 
not uncorrelated. The cross-correlation among various skew-spectra are displayed in Figure [TTI for EPIC and Figure[K)]for Planck. 



CD, 



? (2), 



The dependence of the skew spectra S^' (Ob) and the beam bi (9b) on the smoothing angular scale 8b is being suppressed for brevity. 

The angular bispectrum Bi 1 i 2 i 3 contains all the information at the level of the three-point angular correlation function. These results are generic 
and do not assume any specific form of non-Gaussianity. However, to arrive at specific expressions we will ignore the primordial non-Gaussianity, 
known to be sub-dominant, and concentrate on secondary non-Gaussianity. Ther e is a family of one- point statistics, namely the well-known skewness 
parameters (8b) introduced abo ve, or pseudo-collapsed three-point function dHinsaw et al.ll995h . as well as the equilateral configuration statistics 
dFerreira. Magueiio &Gorskill998h which can all be expressed as linear combinations of the bispectrum terms. The generalized skewness parameters 
introduced above are also all linear combinations of the bispectrum weights but with varying weights. Using one-point statistics has the advantage of 
higher signal to noise but the price we pay is in terms of reduced power to discriminate individual contributions. 

The series expansion for the MFs can be extended beyond the leading order at the level of the bispectrum to the next-to-leading order which 
involves the trispectrum of the temperature field. The lensing induced trispectrum of the CMB will constitute the main next-to-leading order contribu- 
tion. It is also important to realize that measurements of skewness parameters S 1 ' ' (Ob), (8b) and £™ (8b) will not be independent but correlated 
with one another; the level of correlation depends on the noise and beam profile. 

In Figure [2] we have plotted the variance parameters a^(8b) and <r\(8b) for various smoothing beams (assumed Gaussian). The four different 
FWHM that are considered are 8b = 10' , 25' , 50' and 100' respectively. The parameter values only depend on the underlying CMB power spectra and 
the beam as well as the noise. They are used as a normalization parameters while constructing the MFs from the generalized skewness parameters. 
Two different beam and noise levels are considered EPIC (left-panel) and Planck (right-panel). The one-point generalised skewness parameters are 
depicted in Figure-[3]for Planck and Figure-|4]for EPIC. The background cosmology is that of ACDM. 

In Figure |4]the one-point skewness parameters S^'ifib), S^(8b) and S' 2 '(8b) (solid, short-dashed and long dashed lines respectively), are 
plotted as a function of smoothing scale 8b - These parameters are defined in Eq.dl4t. The panels correspond to contributions from different types of 
secondary non-Gaussianity: cross-correlation of lensing and Sunyaev-Zel'dovich effect (left-panel), cross-correlation of ISW and lensing (middle- 
panel); cross-correlation of point source and lensing (right-panel). An experimental set up which is same as EPIC was considered. See Table[T]for 
detail specifications regarding level of noise and beam. The bispectrum used in our calculation is given in Eq.l(5j and the cross-spectra is plotted in 
Figure[3]we plot the one-point skewness parameters for Planck. 



4 THE TRIPLETS OF SKEW-SPECTRA AND LOWEST ORDER CORRECTIONS TO GAUSSIAN MFS 

The skew-spectrum has been studied previously in various cosmological contexts dCooray|[2001al) . e.g. to estimate the bispectrum resulting from 
lensing-SZ correlation. The skew-spectra are cubic statistics constructed by cross-correlating two different fields. One of the fields used is a composite 
field (map) typically a product of two maps either in its original form or constructed by means of applying relevant differential operators. Example 
of such derived maps that we will consider are [Q 2 (Cl)], [VO(J7) • V0(fl)] and [V 2 0(f2)]. The skew-spectra resulting from cross-correlating these 
maps are known as the generalised skew spectra and are related to the three generalised skewness parameters introduced above. At the lowest order, 
the MFs themselves can be constructed using these generalized skewness parameters and contain equivalent information. 

The detection of each individual mode of the primary or secondary bispectrum is still considered challenging. This is primarily due to the 
low signal-to-noise associated with each individual modes. All available information is therefore typically compressed into a single number - the 
skewness. This drastic data compression leads to a significant degradation of the power of the statistic to discriminate between models. 
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Skew-Spectra (ISW x Lensing) 




Figure 6. The skew-spectra sf 0) (0 b ), S*p (O b ) ^ S\ 2) (O b ) for ISWxLensing is plotted as a function of the harmonics I. Notice that skew-spectra corresponding 
to ISW are dominant at smaller is while the ones corresponding to SZ dominate at larger I values. The skew-spectra, being integrated measures, depend on the entire 
harmonic range of the bispectra. The shape of the skew-spectra can play an important role in separating individual contributions of secondary non-Gaussianity. 

The first of the skew-spectra that we will study is the one introduced by dCooravll2001ah and later generalized iMunshi & Heavens! d20ld) . It is 
related to sometimes known as the two-to-one power spectrum and is constructed by cross-correlating the squared map [O 2 (fi)] with the original map 
0(O). The second skew-spectrum is constructed by cross-correlating the squared map [0 2 (f2)] with [V 2 €>(f2)]. Analogously the third skew-spectrum 
represents the cross-spectra that can be constructed using [Ve(fi) ■ Ve(fi)] and [V 2 9(f2)] maps. 

S U^W) S ^ e) ^ S 124^)2m EReal([e]i '" [e2lL) = 127ro^0 ^ Bl hh J ihh b i( e b)b h (e b )b h (0 b ), (18) 

= Tr 2,1 N 2^ ^ (e2 ' V2e) (^) = 1fi 27^ 2^o,L 1 E Real ([ V2Q ]'-[ e2 ]^) 

16-Ka$(6b)<T({0) 16-K(T5{9b)a({6b) 21 + 1 ^ 

= 1G7va 2 a{ g b)a 2 {db) Y, ['(' + + h( - h + ^ + h{ - h + !)] ^1*2 ^M^M^MM^), (19) 

= 87rg 4 (6 , b) ^ [[/(/ + 1) + + 1) - Z 2 (J 2 + l)]ia(/a + 1) + cyc.perm.] Bu^Ju^M^bi^b^eb), (20) 



r _ Ihhh (2l 2 + l)(2l 3 + 1) ( h l 2 l 3 , 



2h + 1 Y (2h + 1)4tt V 

5 W (^) = ^(2« + l)S ; W (e 6 ). (22) 

i 

Each of these spectra probes the same bispectrum Bu 1 i 2 but with different weights. Each triplet of modes specifies a triangle in the harmonic domain 
and the skew-spectra sum over all possible configuration of the bispectrum keeping one of its sides I fixed. 

The extraction of skew-spectra from data is relatively straightforward. The procedure consists of the construction the relevant maps in real space 
either by algebraic or differential operations and then cross-correlating them in multipole space. Issues related to mask and noise will be dealt with 
in later sections. We will show that even in the presence of a mask the computed skew spectra can be inverted to give a unbiased estimate of all-sky 
skew-spectra; presence of noise will only affect the scatter. We have explicitly displayed the experimental beam bi in all our expressions. 

In Figs |5|7| we have presented the three different skew-spectra Si as a function of the harmonics I. The skew-spectra for a generic bispectrum 
is defined in Eq.d 1 81l. Eq .d 19t and in Eq.POt, In Figure[5]we present the skew-spectra corresponding to the SZ effect cross-correlated to lensing. The 
Figure [6] we present the skew-spectra for the ISW effect and Figure [7] shows the skew-spectra for unresolved point sources. The skew spectra are 
sensitive to the beam bi(6t) moreover the skew-spectra at a given I, i.e S^^Ob) depend on the bispectrum Bi 1 i 2 i 3 defined over the entire range of 
modes specified by all possible I values that are being probed. The distinct shape of these individual spectra can be used to study the nature of their 
origin (i.e. primordial or secondary). Specific models of primary non-Gaussianity such as local or equilateral too will have distinct shapes for the 
Sf {Ob) parameters though such contributions will be subdominant for currently allowed levels of primordial non-Gaussianity. 

It is important to stress that these three skew-spectra do not contain completely independent information; the errors associated with them are 
correlated. We next turn to a detailed derivation of signal-to-noise level of these estimators and the level of cross-correlation among these spectra for 
a given observational strategy. The derivations are accurate for near all-sky coverage; for more accurate modeling a computationaly expensive but 
conceptually straightforward Monte-Carlo analysis is required. 
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Figure 7. The skew-spectra (${;), (9j,) and Sp' (Of,) for Point-Source X Lensing is plotted as a function of the harmonics I. The similarity of the underlying 
bispectrum for unresolved point-sources and the SZ effect means the shape of resulting skew-spectra are similar. 



Mission 


e b 






Frequency 


Planck 


7.1' 


2.2 x 10" 6 


0.0349 


143 (GHz) 


EPIC 


5.0' 


8.0 x 10~ 9 


0.002 


150 (GHz) 



Table 1. Parameters used to compute the skew-spectra and the associated scatter for the two different experiments, ongoing Planck I Planck Collaboration 2008) and 
EPIC iBaumann et"aT]|2009h . 



Using the estimator Eq.(|T8jl previous studies have f ocused towards a detec tion of lensing-secondary correlation for individual WMAP frequency 
channels using raw as well as frequency-cleaned maps dCalabrese et~ai]|201f|) . These studies used the KQ75 mask and were limited by the WMAP 
resolution N s id c = 512 and an Z max = 600. No significant evidence for a non-Gaussian signal from the lensing-secondary correlation was found 
in any of the individual bands, 2a and 3a evidence were obtained both for lensing-ISW and lensing-SZ signals in the foreground cleaned Q-band 
maps, respectively. They also found that the point source amplitude at the bispectrum level to be consistent with previous measurements. With higher 
resolution maps available from Planck as well as other future missions such as EPIC it will be possible not only to achieve a cross validation using 
multiple skew spectra, but it should also be possible to reconstruct the topological properties and compare them with the ones obtained in the pixel 
domain. 

A great deal of attention has recently been focused on designing optimal estimators. Indeed this is true that for current generation of experiments 
(WMAP) the mere detection of non-Gaussianity remains a challenging task because of the low ratio of signal to noise. Optimality of an estimator may 
not be a crucial issue for high resolution data from experiments such as Planck, at least for the detection of secondaries, as very high level of signal 
to noise is expected. Attention then will shift to the characterization of non-Gaussianity using an array of estimators to separate different components 
of non-Gaussianity (primordial and secondary) and provide the level of contamination from foregrounds such as point sources. The skew-spectra 
associated with MFs can play a valuable role in this direction. The main advantage of computing the skew-spectra being a direct estimator which can 
be computed using a Pseudo-Cf based approach, and the covariances can be characterized analytically even in the presence of an arbitrary mask and 
non-stationary noise. 



5 ESTIMATOR, SKY COVERAGE AND ERROR ANALYSIS 

The results derived above correspond to a situation in which an all-sky map is available which is free from noise. However, in reality often we have to 
deal with issues that are related to the presence of a mask and (possibly inhomogeneous) noise. Partial sky coverage introduces mode-mode coupling 
in the harmonic domain in such a way that individual masked harmonics become linear combinations of all-sky harmonics. The coefficients for this 
linear transformation depend on specific choice of mask through its own harmo nic coefficients. W e will devise a method that can be used to correct 
fo r this mode-mode cou pling based on the Pseudo-Cf (PCL) method devised by Hivon et al. U2002I) for power spectrum analysis and later developed 



bv lMmshietlirifeOllah for analyzing the skew spectra and the kurt-spectrum jMunshi etalj|201ld) . 

Consider two generic fields A(Q) and B(Cl) and denote their harmonic decompositions in the presence of a mask w(Q) as A\ m and Bi m . Notice 
that the mask is completely general and our results do not depend on any specific symmetry requirements such as the azimuthal symmetry. The fields 
A and B may correspond to any of the fields we have considered above. In a generic situation A and B will denote composite fields and the harmonics 
Ai m and Bi m will correspond to any of the harmonics listed in Eq.ll22b i.e., [0 2 ]; m , [VO • V8]j m and [V 2 0]; m : 
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A lm = / dti YC m (Cl) [w(Cl) A(Cl)]; (23) 



Aim = (— l) m Iihiz ( * ' 2 ' )wi imi Ai 2 , n . 2 . (24) 



liTfli 



Similar expressions holds for B. The above expression relates the masked harmonics denoted by Ai m and Bi m with their all-sky counterparts Ai m 
and Bi m res pectively. In the ir derivation we use the Gaunt integral to express the overlap integrals involving three spherical harmonics in terms of the 
3j symbols dEdmondsll 19681) . The matrix Ii 1 i 1 i 3 encodes the overlap integral defined in Eq.((6}. The expressions also depend on the harmonics of the 
mask wim- If we now denote the (cross) power spectrum constructed from the masked harmonics and denote it by Si(9 b ) and its all-sky counterpart 
by Si(9b) we can write: 



(25) 



sf' B (9 b ) = a^yX^mBim; sf' B (8 b ) =^2M u ,sf' B (e b y, Mu, = ^ £ I 2 ; 

m V VI" 

sf- B (e b ) = ^[M- 1 ]^'^); (ssf' B (e b )ss$' B (e b )) = £ M^(6st' B (9 b )6S A ; B (9 b ))M£,],-, (sf- B (e b )) = s ; A >> 6 ); (26) 

l' LL' 

8sf' B (e b ) = (s A < B {9 b ))~s A < B {9 b y, {a, b} g {e,e 2 ,(ve- ve),v 2 e}. (27) 

The final expressions will be independent of the azimuthal quantum number m due to our assumption of isotropy. In the above derivation we have 
used the orthogonality properties of the 3j symbols. 

It is interesting to note that the convolved power spectrum estimated from the masked sky is a linear combination of all-sky spectra and depends 
only on the power spectrum of the mask used. The linear transform is encoded in the mode-mode coupling matrix Mn, which is constructed from the 
knowledge of the power spectrum of the mask. In certain situations where the sky coverage is very restricted the direct inversion of the mode-mixing 
matrix M may not be possible sue to its singularity and binning may be essential. Based on these results it is possible to define an unbiased estimator 
that we denote by §f' B (9 b ). The noise, due to its Gaussian nature, does not contribute in these estimators which remains unbiased. However, the 
presence of noise is felt in the increase in the scatter or covariance of these estimators (which can be computed analytically): 

[Vt 2 \ = D M-1 ]«' PPk {5v£ 2) SV™) = WlVP^SPPMMZ')' ■ (28) 

I' LL' 

The symbol Sf~' A (9 b ) denotes the power spectrum of the field A(Cl); A(Cl) is a generic field that are used for the construction of generalized skew- 
spectra. The derivation of the covariance depends on a Gaussian approximation i.e. we ignore higher-order non-Gaussianity in the fields. Ci is the 
ordinary CMB power spectra it also includes the effect of instrumental noise and beam Ci(9 b ) — Cfbf(9 b ) + ni. For a survey with homogeneous 
noise, we can write C ; = tlpO-fj where fl p is the pixel area and crjv is the noise r.m.s. In a regime when noise contributions dominate, the MFs can 
be approximated by a Gaussian approximation. The resulting expressions are listed below: 

{6Sf B {9 b )SSf B {0 b )) = ^-j [Sf- A (8 b ) Sf- B (8 b ) + [S A ' B (9 b )] 2 ] ; (29) 

(SSf 1 ' Bl (8 b )SS A2 - B2 (9 b )) = [s Al ' A2 (9 b ) S Bl ' B2 (9 b ) + S Al ' B *(9 b ) S? a ' Bl (0 b j\ ; (30) 

{a u b x ,A2,b 2 }€ {e,e 2 ,(ve- ve),v 2 e}. pi) 



In the next section we will provide detailed explicit expressions for various choices of estimators. 

The two-to-one estimators are from a family of non-Gaussian estimators. Th e three-to-one estimator probes the four-point correlation function or 
equivalently the (angular) trispectrum T/ 1 /^ (L) dCoorav. Li & Melchiorrill2008h . These spectra have been used to probe primordial non-Gaussianity 
beyond lowest order i.e. to separate contributions from tnl and ijnl- The two-to-two spectru m or the power-spectrum of the squared CMB maps was 



found to be useful in the context of studies of weak lensing induced CMB non-Gaussianity (Coorav & Kesden 



120031) . The next order corrections 



the generalized skew-spectra too have been constructed using such an approach dMunshi. Smidt & Cooravll20ld) in the context of studies of MF. The 



to 



Pseudo-C; formalism discussed above is useful for constructing a quadruplets of trispectra. The error estimates for these higher order contributions 
too can be computed using a formalism we outline next. 



6 EXPLICIT EXPRESSIONS FOR COVARIANCES 

As we have already stressed, the estimators we have introduced for the skew-spectra are correlated and do not carry independent information. Their 
correlation structure depends on the experimental beam, noise and sky-coverage. lust as with non-Gaussianiaty, partial sky coverage also introduces 
mode-mode coupling. However using the mode-mixing matrix defined in Eq.J25t it is possible to deconvolve the convolved skew-spectra Si. In this 
section, we list the expressions for the co-variance of various estimators for skew-spectra. The variances or the scatter of the skew-spectra defined in 
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Eq.\2H take the following forms: 

<[< 2 ' W> = {{Sf' e m 2 } - (Sf< e (6 b )f = ^ [Sf- &2 (e b ) Sf-» + [sf-°(6 b )f] , (32) 

([5Sf 2 ^(8 b )f) = ([Sf^(8 b )f) - (Sf^ @ (8 b )f = 0^ [sf' B \6 b ) Sf^(9 b ) + [sf^(e b )f] , (33 ) 

/rrcV 2 e,ve ve//) m2, _ /rc ,v 2 e,ve ve//) u2\ /c <v 2 e,ve ve/ fl >\2 
{[db l (V b )\ ) = (0b) J ) - (6, (0bj) 

/sky r c ,v 2 e,v 2 e// ) n ove-ve.ve ve//) , r a v 2 e,ve-ver a m21 
= 21 + 1 [ 1 { h ' 1 ^ b ' + ' 1 ^ b >> \ ' ( ' 

Here / s k y is the fraction of the sky covered. The expressions for the cross-covariances can also be computed using the similar technique. We express 
the harmonics of the composite fields such as [0 2 ]; m , [VB • VO]; m in terms of harmonics of the [0]; m field using the overlap integral. The final 
equations are derived using Wick's theorem to simplify the resulting expressions. 

« 2 >>b) <5S ; ve ' Ve ' A2e (<M> = ^ [[sf^ e (6 b )][Sf > A2e (<M] + [Sr^ & (9 b )][sf^ & (e b )]\ , (35) 

(6sf-°(8 b ) 5 S f^ e (e b )} = [sf' e2 (6 b ) sr' A2e (e b ) + Sf- e (6 b ) Sf- A2e (e b )} , (36) 

« e ' v9 ' A2e (9 6 ) Ssf' A20 (9 b )) = [Sf 2 ' ve ' ve (0b) sf e ' A2e (6 b ) + sf* A2e (8 b )S7 e - v0 ' A2e (8 b )] . (37) 

The final expressions that we derive are applicable to near all-sky surveys. When a small portion of the sky is covered a sky patch version of our 
calculations can be performed using two dimensional Fourier analysis instead of the spherical harmonic analysis that we use here. Some of the terms 
appearing in these expressions can be expressed in terms of the bispectrum. If we assume that the instrumental noise is Gaussian then there is no 
contribution from noise in these expressions. 

'max 

S ; <e2 ' v2e> (0 6 ) = -^1(1 + 1) Bu lh J ni i 2 bib h b h , (38) 

l i= 2 
i max 

fif ' ve ' ve> (6>b) = B lhh J llll2 [l(l + 1) + h(h + 1) - l 2 (h + l)]bib h b h , (39) 

(;=2 

= - J2 BihhJihhW + !)['(' + 1) + + 1) - h(h + l)]}bib h b h . (40) 



j<ve-ve,A e>^ 



Notice that these expressions are generic, in that they are derived without any specific assumption about the shape of the bispectrum. The rest of the 
terms can be expressed in terms of the power spectrum alone. As is common practice in the literature these results ignore all higher order correlation 
beyond the bispectrum. 



s <e 2 ,ve.ve> ( ^ )= 1 gV h h J V ja, 3 Ji( t + 1) + !) _ + i)]( Cjl ^ + n ll )(C la ^ a +n la ), 

1—9 V / 



(41) 



S 



l i= 2 

2 ^ ( h h I s 



1 21 



l i= 2 
I, 



1 ( o ) fiihiiC'itfi +ni 1 ){Ci 2 bi 2 +n h ), (42) 



g<v-v,v-v) 2 \ ^ ( h h I 



+ tE o ) T hi2iW + V+ Hh + l)-Z 2 (/2 + l)] 2 (C il 6 ! 2 1 +n il )(C i2 6? 2 +n i2 ). (43) 

The remaining terms are scaled input power-spectra: 

5 <e,v»e> = _, ( , + 1){Cibf + ni) . 5 (v 2 e,v 2 e> = f{l + 1)2(C;fe2 + nj); ^e.e = ^ + nj)< (44) 

The derivation of these results follow the same general principle that is outlined in §5\ These expressions are used to compute the cross-correlation 
coefficient among various spectra which are defined below: 



r?(9 b ) = (e b )8S^ (fc)>/V ([SS™ (e b )] 2 )\f{[5Si j) (0b)] 2 ) ; i,j € {0, 1, 2}. (45) 
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Figure 8. The mode-by-mode signal to noise ratio is depicted for the three skew-spectra, (Of,) (left-panel), (Of,) (middle-panel) and (Of,) (right-panel), as 
a function of the harmonics I. Two different experimental set-ups are considered ongoing experiment Planck (solid curves) and EPIC (dashed-curves). The parameters 
defining these experiments are summarized in Table[T] We have assumed an all-sky coverage f s ^ y = 1. The mixed bispectrum that the skew-spectra sample correspond 
to cross-correlation of lensing and thermal SZ effect. All other source of secondary and primordial non-Gaussianity has been ignored. We use Eg. fl29t in conjunction 
with Eg. j31t in our computation of the scatter. 



As before, throughout we have ignored the mode-mode coupling. The coefficients of cross-correlation ty, are independent of the sky-coverage / s ky 
The signal to noise for individual modes for a given spectrum on the other hand can be expressed as: 



[S/tyfm = \J{[s i ; ) (8 b W)/{[ss ( ;\e b )]^) i e {o, 1,2}. (46) 

The cumulative [S/N] = X^!=2 X [S/N]^ is tabulated for individual experiments in Table-[2]for Planck and EPIC. 

We have also computed the signal to noise ratio for individual modes using these expressions for various skew spectra. These results are plotted 
in Figure|9]for EPIC as well as for Planck in Figure ([8). The cumulative signal-to-noise as expected is higher for EPIC due to higher sensitivity. The 
signal to noise for ISW decreases sharply at higher I and peak at lower / on the other hand the signal to noise for SZ and unresolved point sources 
peak at a much higher angular frequency. Among the three skew-spectra we have considered, the skew spectra • (0b) was found to have higher 
signal to noise compared to (Of,) and sjy (8b)- While the lowest order skew-spectra 5^°' (Of,) is dominated mostly by cosmic variance the other 
skew-spectra, SP (0b) is maximally affected by the instrumental noise. The information content is not independent for the different skew-spectra; 

their cross-correlation coefficient provides a succinct measure of this lack of dependence. 

To correct for the effect of a mask and the noise we will follow the Pseudo-Cf (PCL) method devised by lHivon et a l. (2002) for power spectrum 



analysis and later developed by iMunshi et al.l d2011al) for analyzing the skew spectra and the kurt-spectrum fMunshi et al] 201 id) . and the cross- 
correlation coefficient provides a valuable indicator of their independence. 

The signal to noise of estimates of one-point generalised skewness parameters = X^=2 X (^ + l)^t is given by 

[E!=2 X ( 2l + l) 2 {[5S (i) /S (i) } 2 )^ . The corresponding numbers for Planck and EPIC are presented in Tableg] 

It is possible to introduce a filter function wi 1 i 2 in the definition of the skew-spectra. Choices include, sharp cutoff in the I space to avoid the 
affect of noise at high I, to optimal filters that maximizes the signal to noise for a given resolution i m ax- Clearly such option will invariably improve 
the statistical significance. The filter functions can be of further interest if the bispectrum is more pronounced for certain triangular configurations. 
Another potentially useful application of a filter function is to filter out a specific configuration. However in such case the skew-spectra will not have 
any direct relation with the topological properties of the original map. 



7 CONCLUSION AND DISCUSSION 

Non-Gaussianity is in itself a poorly defined concept, in that there is no unique approach that can be adopted to describe or parametrize an arbirary 
form of non-Gaussianity in a complete manner. In order to quantify non-Gaussianity as fully as possible it is therefore essential to deploy a battery 
of complementary approaches each of which exploits different statistical characteristics. Each such technique will have a unique response to real 
world issues such as the sky-coverage (observational mask), beam and instrumental noise. Any robust detection therefore will have to involve 
a simultaneous cross-validation of results obtained from independent methods. The most common characterizations of non-Gaussianity involve 
studying the bispectrum, which represents the lowest-order departure from Gaussianity; higher order non-Gaussianity can be studied using its higher 
order analogues i.e. the multi-spectra. 

By contrast the topological estimators (MFs) that we have studied here carry information to all orders, though in a collapsed (one-point) form. 
Analytical results for MFs for a Gaussian field are well understood, and form the basis of non-Gaussianity studies jTomitJl986h . There have been 
several previous studies on extraction of the MFs from the CMB data that rely either on simplification of radiative transfer using the Sachs-Wolfe 
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EPIC (signal to noise) 




Figure 9. The signal to noise ratio for the three skew-spectra, (0^) (left-panel), Soffit,) (middle-panel) and (6^) (right-panel), are depicted as a function 
of the harmonics I. Two different experimental set-ups are considered ongoing experiment Planck (solid curves) and EPIC (dashed-curves). The parameters defining 
these experiments are summarized in Table [T] We have assumed an all-sky coverage f s ^ y = 1. The mixed bispectrum that the skew-spectra sample correspond to 
cross-correlation of lensing and thermal SZ effect. All other source of secondary and primordial non-Gaussianity have been ignored. 



(Planck,EPIC) 


SZ 


ISW 


PS 


(S/N) 


(5.0, 1137.4) 


(1.0,216.0) 


(0.5,209.0) 


(S/N) 


(24.0, 1354.9) 


(62.2,420.3) 


(4.3,552.0) 


(S/N) 


(19.7, 1328.8) 


(31.8,246.5) 


(1.7,421.0) 



Table 2. The cumulative signal to noise (S/N) for Planck and EPIC surveys, a re shown for the three one-po int skew-s pectra. Parameters use d to compute the skew-spectra 
and the associated scatter for the two different experiments, ongoing Planck i Planck Collaboration 2008) and EPIC (Baumann et al. 2009). 



limit faikage et al.l2006h or using a perturbative approach based on a series expansion of the MFs that can be studied order by order (Matsubara I2003L 
Il994l) . The MFs have also been studied using elaborate computer-intensive non-Gaussian simulations dKomatsu et al ■l2003l ; ISpergel et al.l2007l) . Most 
of these studies were done using a specific model of non-Gaussianity, namely the local model of primordial non-Gaussianity which is parametrized 
by the well known parameter /nl- 

The main motivation behind the present study has been to to extend such methods to secondary non-Gaussianities which have not been studied 
before in the context of morphological statistics analytically. The increase in sensitivity of CMB experiments and near all-sky coverage along with 
wide frequency range means the study of non-Gaussianities will be feasible in the very near future. Moreover, in the currently favoured adiabatic 
CDM models it is expected that the contribution from primordial non-Gaussianity is negligible and the main contribution to non-Gaussianities comes 
from secondaries. The secondary no n-Gaussianity signal are ass o ciated with large sca le structure contributions and through various mode coupling 
effects such as gravitational lensing l lGoldberg & Spergelll999al lbl lCoorav & HukOOd) . Our primary aim in this work has been to study how well we 
can probe the secondary signals from mode coupling using morphological descriptors. 

One of the main difficulties faced by one-point estimators S^'(6b), which also affects the MF-based estimators 14 (V), is their inability to 
differentiate among various sources of non-Gaussianity. The triplets of skew spectra £Tj (9b) that we have introduced can be used to separate out 
contributions from various s econdaries as well as to pro be and constrain any foreground residuals left from the component separation step of the data 
analysis chain. Generalizing iMunshi & Heavens! J2010h we have introduce a set of triplets of skew-spectra which can be extracted from any realistic 
data. These skew-spectra do not compress the available information from a bispectrum to a single number, and their shape can help to distinguish 
among various sources of non-Gaussianity. Exploiting the perturbative expansion of the MFs, we showed that at the leading order of non-Gaussianity 
the MFs are completely specified by the knowledge of the bispectrum. Our results are most naturally defined in the harmonic domain. Comparison of 
MFs extracted using harmonic approach can be cross-compared with more traditional approach in the real space as an useful consistency check. 

The methods based on the skew-spectra that we have presented are simple to implement once the derivative fields [VO ■ VO] o r [V 2 €>] are 
const ructed. We have shown that this can be implemented in a model-independent way. Our method is based on a Pseudo-C; approach dHivon et al.l 
l2002h and can handle arbitrary sky coverage and inhomogeneous noise distributions. The Pseudo-C^ approach is well understood in the context of 
power spectrum studies, and its variance or scatter can be computed analytically. We provide generic analytical results for the computation of scatter 
around individual estimates. We also provide detailed predictions on how they are cross-correlated. In our method, it is possible indeed to go beyond 
the lowest level in non-Gaussianity to include the contribution from trispectrum. The main contributions in frequency-cleaned CMB maps will be 
from lensing of the CMB, though it is expected that such corrections will be sub dominant at least in the context of CMB data analysis. 

We conclude by pointing out that the MFs do not probe the full bispectrum, but involve only weighted sums of modes and are thus equivalent to 
the three generalised skewness parameters we have used. We have also defined three generalized skew-spectra associated with each of these skewness 
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(Planck,EPIC) 


sz 


ISW 


PS 


(S/N) 


(3.8,503.2) 


(0.3,39.3) 


(0.4, 53.4) 


(S/N) 


(5.8,1299.0) 


(6.4 x 10" 2 ,70.9) 


(0.6,529.1) 


(S/N) 


(1.2,625.8) 


(1.4 x 10~ 2 ,21.1) 


(0.1, 178.2) 



Table 3. The cumulative signal to noise (S/N) for Planck and EPIC surveys for the one-point cumulants SW (57), defined in Eq.{46}, are shown for each sk ew-spectra. 
Parameters used to co mpute the skew-spectra and the associated scatter for the two different experiments, ongoing Planck IPlanck Collaboration! l200j) and EPIC 
iBaumann et al.l2009l) . We have assumed / s k y = 1 in our calculations. 



Planck (Cross - correlation) 




Figure 10. Same as the previous Figure but for Planck noise. The expression for scatter Eq. d32t -Eq. )34t and cross-correlation Eq. )35t -Eq. )37t has two contributions. In 
each of these expressions there are terms which depend on the bispectrum and there are terms which can constructed from power spectrum alone. For Planck noise we 
found that the expressions for the scatter as well as the cross-correlation are entirely dominated by the terms which depend only on the power spectrum thus making the 
coefficient ry independent of the type of underlying bispectrum. 



parameters. In this sense, the study of these skew-spectra can replace the study of MFs. The skew-spectra we have introduced can all be probed for 
arbitrary mask and noise. Ubiased estimators can also be constructed which can work in the presence of partial sky coverage and inhomogeneous 
noise. Their variance can also be computed analytically; thereby avoiding the use of non-Gaussian Monte-Carlo simulations completely. Finally, the 
MFs can be constructed from the knowledge of generalized skew-spectra and can be compared with the results from real space analysis. The triplets 
of generalized skew-spectra can be used to separate individual components of NGs using their shape information. From our analytical results of 
cross-correlation, we find that in the absence of noise, e.g. experiments such as EPIC, the skew-spectra are highly correlated, more so for higher I 



EPIC (Cross - correlation) 




500 1000 
I 




500 1000 
I 



0.5 



-0.5 



PS 



500 1000 
I 



Figure 11. The information content of the skew-spectra are not completely independent. The level of cross-correlation among various estimator is encoded in the 
coefficient of cross-correlations ry, defined in Eq. )45t . The cross-correlation coefficient roi, ro2 and ri2 are plotted as a function of harmonics I. The noise level 
correspond to that of EPIC (see Table[T}. The correlation structure refelects the underlying spectra 0i as well as the level of noise. The cross-correlation is similar for 
SZ and PS. 
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values. The correlation coefficients are typically in the range r = 0.5 — 1 for a Planck type experiment. The cumulative signal-to-noise ratio, in a 
Planck type experiment, for bispectrum corresponding to the ISW and SZ and lensing cross-correlation reaches O(10). An improvement of about 
two orders of magnitude can be expected with experiments such as EPIC. 

Throughout we have ignored the presence of primordial non-Gaussianity which is expected to be subdominant. Nevertheless, it can be incorpo- 
rated. Individual skew-spectra from different underlying bispectrum can essentially be combined to construct the total skew-spectra which means that 
our results can straightfowardly be generalised to incorporate specific models of primordial non-Gaussianity. 

The generic results derived here ar e also applicable to othe r areas in cosmology a nd have indeed been explored recently. Examples include 
the analysis of galact ic redshift surveys l lPratten & Munshill2012t) , weak lensing surveys dMunshi et alj|2012ah and the frequency cleaned SZ maps 
l lMunshi et al.ll2012bT) . The results presented here can be extended beyond the analysis of temperature maps, e.g. to analyze polarisation maps, by 
extending the spin-0 calculations to spin-2. Such results can furnish useful probes for tje characterization of morphology of reionization in three 
dimensions. 

A few comments are in order about the comparison of our estimators with the so-called optimal estimators. The motivation to construct an 
optimal estimator is to improve the signal-to-noise of detection which is important in case of weak signals such as the primordial non-Gaussianity. 
The main motivation in this paper has been to reconstruct the topological properties of the CMB map going beyond Gaussianity, in the harmonic 
domain; in particular due to the contributions from secondary lensing cross-correlation which will be detected with high signal to noise with the 
proposed CMB surveys such as EPIC. 

In addition to the primary and secon dary non-Gaussianity, cosmic defects such as textures or cosmic strings jAlbrecht. Battev & Robinson! 19991 ; 
ICruz etal.l2007l ; |Regan & Shellard201(j|) also leave non-Gaussian footprints in CMB maps which can be detected by the change in topological nature 
of the maps. The estimators we have presented here may have relevance in such investigations. A detailed study will be presented elsewhere. 

At the level of the bispectrum the effect of lensing can only be studied through its cross-correlation with other secondaries. However weak 
lensing is also independently responsible for a next order correction to MFs through its effect on the trispectrum; the signal-to-noise is expected to be 
low. 

The signal-to-noise of the skew-spectra fo r secondary-lensing cross-correlation bispectrum is comparable to that of the skew-spectra of 
frequency-cleaned SZ maps (Munshi e t alj|2012bh . However the secondary skew-spectra are much higher compared to skew-spectra associated with 
primary skew-spectra unless we assume a rather high value for the /nl- 



8 ACKNOWLEDGEMENTS 

DM and PC acknowledges support from STFC standard grant ST/G002231/1 at the School of Physics and Astronomy at Cardiff University where 
this work was completed. DM would like to thank Joseph Smidt, Geraint Pratten, Asantha Cooray, Shahab Joudaki and Erminia Calabrese for very 
useful discussions. 



REFERENCES 

Acquaviva V., Bartolo N., Matarrese S., Riotto A., 2003, Nucl. Phys. B667, 1 19 

Adler R. J., 1981, The Geometry of Random Fields, Chichester: Wiley 

Albrecht A., Battye R.A., Robinson J., PRD, 1999, 59, 023508 

Babich D., 2005, Phys. Rev. D72, 043003 

Bartolo N., Matarrese S., Riotto A., 2006, JCAP, 06, 024 

Bartolo N., Komatsu E., Matarrese S., Riotto A., 2004, Phys.Rept. 402, 103 

Baumann D. et al. [CMBPol Study Team Collaboration], AIP Conf. Proc.2009, 1 141, 10 arXiv:081 1.3919 

Birkinshaw M„ Phys.Rept, 1999, 310, 97 

Bock J. et al. 2008, arXiv:0805.4207 

Bock J. et al. 2009, arXiv:0906.1 188 

Boughn S., Crittenden R., 2004, Nature, 427, 45 

Calabrese E., Smidt J., Amblard A., Cooray A., Melchiorri A., Serra P., Heavens A., Munshi D. 2010, PRD, 81, 3529 

Canavezes A., et al., 1998, MNRAS, 297, 777 

Castro P. G. Phys.Rev. 2004, D67, 044039 

Chen X., Advances in Astronomy, 2010, 2010:638979 

Coles P., 1988, MNRAS, 234, 509 

Cooray A.R., Hu W., 2000, ApJ, 534, 533 

Cooray A., 2001a, PRD, 64, 043516 

Cooray A., 2001b, PRD, 64, 063514 

Cooray A , Sheth R., 2002, Phys.Rept., 372, 1, 

Cooray A., 2002, PRD, 65, 103510 

Cooray A., Kesden M., 2003, New Astron., 8, 231 

© 0000 RAS, MNRAS 000, 000-000 



16 Munshi, Coles & Heavens 



Cooray A., 2006, PRL, 97, 261301 

Cooray A., Li C, Melchiorri A., 2008, PRD, 77,103506 

Copi C, Huterer D., Schwarz D., Starkman G., 2007, PRD, 75, 023507 

The COrE Collaboration, 201 larXivl 102.2181 

Cruz, M., Turok N., Vielva P., Martnez-Gonzlez E., Hobson M, Science 318 (5856) 1612 

Edmonds, A.R., Angular Momentum in Quantum Mechanics, 2nd ed. rev. printing. Princeton, NJ:Princeton University Press, 1968. 

Falk X, Madden R., Olive K.A., Srednicki M., 1993, Phys. Lett. B318, 354 

Ferreira P.G, Magueijo J. & Gorski K.M., ApJ, 1998, 503, 1 

Gangui A., Lucchin E, Matarrese S„ Mollerach S., 1994, ApJ, 430, 447 

Gleser L., Nusser A., Ciardi B., Desjacques V., 2006, MNRAS, 370, 1329 

Goldberg D.M., Spergel D.N., 1999a, PRD, 59, 103001 

Goldberg D.M., Spergel D.N., 1999b, PRD, 59, 103002 

Gott J. R„ Mellot A. L., Dickinson M., 1986, ApJ., 306, 341 

Gott J. R., et al., 1989, ApJ., 340, 625 

Gott J. R., et al., 1990, ApJ., 352 

Gott J. R., Mao S., Park C, Lahav 0., 1992, ApJ., 385, 26 

Hadwiger H., 1959, Normale Koper im Euclidschen raum und ihre topologischen and metrischen Eigenschaften, Math Z., 71, 124 

Hanson F.K., Lewis A., 2009, PhRvD, 80, 063004 

Hikage C, et al., 2002, Publ. Astron. Soc. Jap., 54, 707 

Hikage C, Komatsu E., Matsubara T., 2006, ApJ., 653, 1 1 

Hikage C, et al., 2008, MNRAS, 389, 1439 

Hikage C, et al., MNRAS, 2008, 385, 1613 

Hikage C, Komatsu E., Matsubara T., 2006, Astrophys. J., 653, 1 1 

Hinsaw G., Banday A.J., Bennett C.L., Gorski K.M., Kogut A., 1995, ApJ, 446, 67 

Hivon E., Gorski K. M., Netterfield C. B., Crill B. P., Prunet S., Hansen E, 2002, ApJ, 567, 2 

Hoftuft J., Eriksen H.K., Banday A.J., Gorski K.M., Hansen F.K., Lilje P.B., 2009, ApJ, 699, 985 

Hu W., Scott D., Silk J., 1994, PRD, 49, 648 

Joudaki S., Smidt J., Amblard A., Cooray A., 2010, JCAP, 08(2010)027 

Kaminkowski M., Spergel D., 1994, ApJ, 432, 7 

Kamionkowski M., Smith T.L., Heavens A., 2011, PRD, 83, 023007 

Kerscher M., et al., 2001, A&A., 373, 1 

Komatsu E., et al., 2003, ApJS, 148, 119 

Kofman LA., Starobinsky, 1985, Sov. Astro., Let., 11, 271 

Larson D. et al, ApJS, 2011, 192, 16 

Luminet J.P., Weeks J., Riazuelo A., Lehoucq R., Uzan J.P., Nature, 2003, 425, 593 
Maldacena J.M., 2003, JHEP, 05, 013 

Martinez-Gonzalez E., Sanz J.-L, Silk J., 1990, ApJ, 355, L5 

Matsubara T., 1994, ApJ, 434, L43 

Matsubara T., 1995, ApJS., 101, 1 

Matsubara T., Jain B., 2001, ApJ, 552, L89. 

Matsubara T., 2002, arXiv:astro-ph/0006269v2 

Matsubara T., 2003, ApJ, 584, 1 

Matsubara T., 2010, PRD, D81, 083505 

Matsubara T., Yokohama J., 1996, ApJ, 463 

Matsubara T., 2010, PRD, D81, 083505 

Matsubara T., Jain B„ 2001, ApJ, 552, L89 

Melott A. L., 1990, Phys. Rep., 193, 1 

Moore B., et al., 1992, MNRAS, 256, 477 

Mollerach S., Gangui A., Lucchin E, Mataresse S., 1995, ApJ, 453, 1 
Mukhanov V.F., Feldman H.A., Bandenberger R.H., 1992, PRD, 215, 203 
Munshi D., Souradeep T., Starobinsky, ApJ, 1995, 454, 552 
Munshi D., Melott A.L., Coles P., 2000, MNRAS, 311, 149. 
Munshi D., Heavens A. 2010, MNRAS, 401, 2406 
Munshi D., Smidt J., Cooray A., arXiv: 101 1.5224 

Munshi D., Heavens A., Cooray A., Smidt J., Coles P., Sena P., 201 1, MNRAS, 412, 1993 
Munshi D., Valageas P., Cooray A., Heavens A., 201 la, MNRAS, 414, 3173 
Munshi D., Coles P., Cooray A., Heavens A., Smidt J., 201 lb, MNRAS, 410, 1295 



© 0000 RAS, MNRAS 000, 000-000 



Secondary Anisotropics and Minkowski Functionals 17 

Munshi D., Smidt J., Heavens A., Coles P., Cooray A., 201 lc, MNRAS, 411, 2241 
Munshi D., Etching T., Heavens A., Coles P., MNRAS, 201 Id, 416, 629 
Munshi D., van Waerbeke L., Smidt J., Coles P., MNRAS, 2012a, 419, 536 
Munshi D., Smidt J., Joudaki S., Coles P., MNRAS, 2012b, 419, 138 
Natoli P., et al., arXiv:0905.4301 

Novikov D., Schmalzing J., Mukhanov V. R, 2000, Astron. Astrophys., 364 

Park C, et al., 2005, Astrophys. J., 633, 1 1 

Planck Collaboration, arXiv:081 1.3919 

Pratten G., Munshi D., 2012, MNRAS, 423, 3209 

Regan, D. M., Shellard, E. P. S., 2010, PRD, 82, 063527 

Roukemia B.F., Lew B., Cechowska M., Marecki A., Bajtlik S., 2004, A&A, 423, 821 

Sahni V., Sathyaprakash B.S., Shandarin S.F., ApJ., 1998, 495, L5. 

Sato J., etal., 2001, ApJ, 421, 1 

Salopek D. S., Bond J. R., 1990, PhRvD, 42, 3936 

Salopek D. S., Bond J. R., 1991, PhRvD, 43, 1005 

Schmalzing J., Gorski K. M., 1998, MNRAS, 297, 355 

Schmalzing J., Diaferio A., 2000, MNRAS, 312 

Smidt J., Joudaki S„ Serra P., Amblard A., Cooray A., 2010, PhRvD, 81, 123528 

Smith K.M., Zahn O., Dore O., 2007, Phys. Rev. D, 76, 043510 

Spergel D.N. et al., 2007, ApJS, 170, 377 

Szapudi I., Szalay A.S., ApJ., 1999, 515, L43 

Taruya A., et al., 2002, ApJ, 571, 638 

Tomita H., 1986, Progr. Theor. Phys., 76, 952 

Verde L., Spergel D.N., PRD, 65, 043007 



APPENDIX A: SPHERICAL HARMONICS 

The completeness relationship for the spherical harmonics Yj m (f2) is given by: 

Y / Y lm (fi)Y l * rn (n') = s 2D (Ci-£i'). (Ai) 

Im 

The orthogonality relationship is as follows: 

dCl F im (n)F; m , (fi) = <$><5* m ' ■ (A2) 
Here 820 and 8 K represents the Dirac's 2-dimensional delta function and Kroneker's Delta function resepctively. 

APPENDIX B: 3J SYMBOLS 

The following properties of 3j symbols were used to simplify various expressions. 

£( 2I ' + 1 )( ^ m 2 m s )(mi m' 2 L ) = 6 ^< 6 **<> (Bl) 



/ \ 8 K 8 K 

/ h h h \ I h h I3 \ _ hi' 3 m 3 m' 3 _ 
I mi ni2 m3 J \ mi m 2 m' 3 J 21^ + 1 ' 



m]ra2 



I I I' \ (-1)' 



^ ( m -m ) = ^^'o; (B3) 
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